function f=find_fhat_w(w,W_i,h_w,k_type)

n=size(W_i,1);
f=NaN*w;

for jj=1:size(w(:),1)
	K_w=kfunc((W_i-ones(n,1)*w(jj))/h_w,k_type);
	f(jj) = (1/(n*h_w))* sum(K_w,1);
end
	%
